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j^bstract-The optic disk is Iheeiilpoint of the retinal nerve fibers 
from the eye, and the entrance and exit point for retinal blood 
vessels. It is a brighter region than the rest of the ocular fundus 
and its shape is approximately round. The location of the optic 
disk is crucial in retinal image analysis, for example, as a 
reference to measure distances and identify anatomical parts in 
retinal images (e.g. fovea), for blood vessel tracking, and many 
others. Previous work on the optic disk has mainly focused on 
locating its center only. The aim of this study is accurately 
locating the optic disk and well defining its contour fully 
automatically. The proposed algoiithm is based on 
morphological operations and snakes for precise optic 
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snake must be close to the 
eiges to the wiong place, 
suitable initial snake is a 



I.&tp,odu:tio^] 

R^nal image is a hi^y specialized form of nodical 
imaging. A fundus carr»:a is a spedali^id \ov/ powa" 
rrridoscope with an attached caniE^B. Its optical desgn is 
baaed on the indirect ofJithalnioscope. The photografJier can 
visualize the back of the eye by focusing light through the 
comea, pi^sl and lens. 

Most of the trials on the optic disk mainly focused on 
locating its caita" only [1, 2, 3]. Cox et al. [4] used gr^ levd 
information around Ihe vicinity to automatically extract the 
boundary of the optic di^ with initial apprmmate location 
givKi 1^ the user irput Osar^ et al. [5] build Iheir method 
on terrq:Jate matching and cfynamic contour to deteimine the 
boundary of the optic di^ 

The proposed algorithm consists of different image 
processing techniques starting with the digital color retinal 
ima^ and the final result is fully automated piedse 
localization of the optic didi boundary. 

n.NfeTHODOLOGY 

Retina! Image Data Acquisition 

The c^tured fundus images were obtained using a 
TOPCON TRC-501A non-rnydiatLic fundus cameg, with a 
Soriy CCD video carrera attadied to system Tlie captured 
fundus ima^isa570x57Clx24 bits colordigital image. Fig. 1, 
shows an exanple. 

RGB and HIS CohrModels 

The RGB color model (where R, G, and B are abbreviaied 
from the colors Red, Green, and Blue respectively) is used, in 
this woit, to display ihe color retinal image, and in ihe 
background elimination step, while the HSI color model 
(where H, S, and I are abfareviations for Hue, Saturation, and 



Intensity respectively) is used in the all subsequent processing 
stages. The RGB and HSI have an invertihie relation between 
them [6], 

Background Elimtnotion 

The optimum threshold technique [6] is applied to the red 
component image, to make iJie i \e\s belongng to the f ndus 
image region and the pntels helo ging to the dart surrou d j 
region separable. The ophn u thresl old s ca! ulated us j 
the following iterative aljonthn 

1. Assuming no k owledge abo t tie exa t locaho of 
fundus region, cons de as a f "st approx n ation that the 
four comers of liie mage co ta la kgro nd p xels o Ij 
and the ren a de nta s f m d s prxels 

2. At ste( t o p te / g and ^ as the eai backjro md 
and fi nd s jra le el "espech ely whe -b segn entation 
into ba k J "o u 1 a 1 f m lus at step t s def ed by the 
thresh Id alue T detem ned n Uie p evious step 



if^background _ pixels 

#junaus_ pi\as 
wbiece: lR(x,y) is the red intensity of the pixel (x.y) in the 
color image. 

Set T = w^ch provides an i^xiated 

background/fundus distinction. 

If T -T < epson thai lialt otherwise go to step 



(2). WhaB qDson is a small numba" like oi 






Locd CoTtra^Ertiancerreit 

Because ihe fundus is not a plane (spherical face), 
illumination is uneven. The brightness of objects i 
background in central region is usually higher than that in 
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s TO nd J "ej on A.s a "es t t e o hast of he f d s 
n ages te Is h> d mimsh as the d s ance of a p xel fro the 
e te of the aje creases Hence a prep Dcess g co tras 
enl ancement s neede 1 to pdu e th s effect 

Let/ yj be the nte sly mage of the olo "etjnal n aje 
Consder a s b mage w n low w(\]!) e tered o a | el 
located at ( y) Suppose tl at /„, a 1 / are the n ax 
and n n niensties of the wl ole n age Denote he 

meai a 1 s andand dev atio of the tens ty w 1 Why 
\i (f) ai d (y (f) respectively The adaj h e contras 
e 1 an eme t tra sfo aho [ ] s def ne 1 1 y 

'.(f)-i>„(f„.)] 



/(x,j)^g(x,y)-255 

wiieire the agmoid function is 



[(>.(/„)-(>.(/„)] 



it) 1+exp ' 
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After ai ply j tl s transfo mat on tl £ 
sn all te s ty standard le alio s (j oo 
more detals However tl s w 11 also 
Hence ^D Ga ss an s oo 1 g f te [6] 1 as 1 ee ; 
onder to ed ce the no se 1 efore the local 

CO \o tio ask 1 at ap o a 
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Region of Intel est (ROD Identtfitotion 

The ophc disk is tjpically found in a region that occupies 
approximately 80x80 pixels The appearance is characterized 
by a relatively rapid vanahon in intensity because the "dait" 
blood vessels are beside Ihe "bright" nerve fibers. The 
adjacent pixels intensity vanance is used for initial detection 
of Ihe optic disk center [2] 

Let (jfx y) be the resultant intensity image of previous 
processing step (contrast enhancement) and w(x,yj be a sub 
image window centered on a pixel located at (x,y). A variance 
image is formed by the transformation: 

Pf^.S)-MjQ^)-MliS) (5) 

An image q(x,yj of the average variance within the sub 
images is then oblained as: 

The point (x^,i/i) with the maximum of this image is taken as 
an initial estimation of the center of the optic disk. Then a 
square region of interest (ROI) is constructed such that ihe 
point iXi,i/i) is in its center The size of this sub image is taken 
110x110 pixels to ensure containing Ihe entire area of the 
optic disk. 
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a Est2 nation of e Uph D kCo o 

Tl e Hough tra sfo [6] s a te h que wl cl an be se 1 
to sola e features of a partic lar si a| e w h an n age 
Be a se req "es tl at li e les "ed featu bs be spec f ed n 
some paran etnc fo m the c re lar Hough transfo m can 1 e 
used f o tl e tial approx ate letection of tl e o| ti 1 sk as 
a neariy c rulai olject n the fu lus aje Tl e n otivat g 
dea 1 ehind 1 e Houjh te h que fo c rcle tl e 
ai p "o mat on of Ihe oj tic d sk lete tio s that each put 

eas re ent (e j cooid ate | o t) d ates ts contnbution 
to a glol all cons ste t solution A c rcle can 1 e des nl ed 1 
the paiametnt, equation. 

(x-xtft(y-s,f.r'' (7) 

v\terB \i and y^ are the coordinates of the center of the circle 
ad r s li e ^dius. The transform is implemented by 
q ai hzing tl e Ho igh parameters space into finite intervals or 
ace ulator cells (i.e. a multidimensional array). As Ihe 
algontl m ns each (x,,ij) in the estimated "ROI" is 
transfo "med nto a discretized fxj , yj ,rj and the accumulator 
cells wl cl sat sfy Eq. (7) are incremented Peaks m the 
ace ulator array represent strong evidence that a 
coiiespo di J c T;le exists in ihe image. Computational load 
an 1 e red ce 1 by using the edge's gradient direction vector 
Tl e ce ter of a c rcle lies on a line along the gradient vector 
d e t on Tl s reduces the update to two points at each radius 
[ B] 

How eve tl e actual optic disk contour is not pure circular 
in shape. So, to identify accurately tiie actual contour of Ihe 
optic disk, active contour models or snakes technique is 



Blood Vessels Removal 

Active contour mcxtels are often called aiakes because Ih^ 
^^Dear to dithes" across images (a phKiomraion known as 
l^'stereas). This rr^hod works on a gradiait image and lock 
onto horrogeneous regions enclosed by strong gradieait 
information. In the fundus image and to use aiakes for 
accuratdy locating the optic didt contour, this task is made 
exlien^y difficult sinoe Ihe optic di^ region is invariatJy 
fragnented into mullipJe regions by Hood vessds This is 
vAjy furtiiea" processing is required before ^^lying Ihe snake 
mdhod. 

The 'TIOI" is processed I5' two niDrphological operations 
to remove tiie regions of Hood vessels that are found on Ihe 
optic di^ to create fairiy constant region. This is done by 
^^lying morphological closing opoiation, i.a a dilation 
[xocess [6] to first remave the blood vessds aiKl then an 
aDaon process [6] to restore the boundaries to liieir former 

P e seLo a! latio o/tJeOpScD k 

Active CO tou etl od [^ 10] hin s the piable of edge 
detection a 1 I k j to an ene"gy n m zaho piablem 
Typically the tra litional way of e Ige detecho conssts of 
ai plj ng a g "ad e t o| erator to Ihe mage and tl e 1 k j tl e 
fou 1 often no conb ous edjes The act vecontou o lei 
works the other way around. First a continuous initial contour 
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or curv e model is fonned (the urcle of optiL disk 
appioximahon found previouslv bj the Hough fransform 
method) Then interna! and external forces act on the model 
drawing the curve to an equilibnum position which will be 
the accurate ophc disk contour 

Because of the way the active contour models beha\ e 
while minimizing their energy they are also called snakes 
Fig 3 shows a snake with its ends joined so that it fonns a 
closed loop. Over a senes of time steps the snake moves into 
alignment with the nearest salient feature (in this case an 
edge). The contour is influenced by the following forces: 

1. Internal forces: internal constraints give Ihe model tension 
and stiffness. 

2. External forces: external constraints come from high-level 
sources such as human operators or automatic 
initial izalion procedures. 

3. Image forces: image energy is used to drive Ihe model 
towards salient features such as light and dark regions, 
edges, and terminations. 

The total energy of ihe model (Es„a(a) is given by the sum of 
the eneigy for the individual snake elements (vertices 'v'): 

E..... - Je,«™i (») + Je„.™i (»)+ Je,„.„ (») (a 

The internal ena:gy of a aiake danent is defined as: 

B«,r« (V).2(s)|v,(sf + /;(s)|v„(s)f (9) 

This energy contains a first-order term controlled by a(s), and 
a second order term controlled by ^s). The first-onder term 
makes the snake contract like an elastic band by introducing 
tension, while the second order term makes it resist bending 
by produang stiffness Adjusting the weights ais) and p(s} 
controls the relative importance of the tension and stiffness 
terms 

The external energy term is used to control attraction or 
repulsion forces iJiat drive active contour to or from specified 
features These forces are easily found by differontiabon. The 
gradient vector flow [9,10] provides the externa! force 
npcessary to push the initial circular active contour toward ihe 
actual [ ontoui of the optic disk. 

III. Results & Discussion 

Itisfound that the red conponent of the color retinal inage 
is tlie mast siiitaHe for the background dirnination st^. Fig. 

4. is the resultant lanary image after the background 




n 

Fig. 4, The resiilt of the background elimination step. 

dimination of the color fundus image dwwn in Hg. 1. 

Aft«: dirrinating the Hack background fiDm the fimdus 
image, the color retind image is convated to an intaiaiy 
image {'T' conponent of the HSI rr»dd). Thai, the contrast 
of Ihe intensity image is aihanced by ^fJying Ihe locdly 
additive contrast aihancerrfint tian^ormation rreihod. The 
intaiaty images bdbre and after contrast eiihaiK:emKit are 
diown in Hg. 5. 

As diown in Hg. 5 (b), the optic didi is in a region of 
neariy 8QxB0 pixds and its appearance is characterized by a 
rdativdy r^id variation in intaialy because the "dark" Hood 
vessds are beade Ihe "bri^t" nerve fibos. Uang this 
fJienomenon and Equations (5, 6) the 'HOI" is localized from 
fXDcesang the image shown in Hg. 5 (b). Hg. 6, dioMS the 
result of the localization of "ROI" (labded by fcJack square) 
sup«iiiposed on the original color retinal image of Hg. 1. 

The rest of the procesang steps are ajplied only to the 
resultant 'TIOI" to reduce the fxocessing tirre and not to 
genante a false initial contour for the snake^ espeddly in the 
images of paliaits suEfaing from "diabdic retinopathy". 

The step of ddiamining a circular ^prcodmation of the 
optic disk using the drcular Haugh transform is pecforrrKd to 
obtain an initial contour for the aiake rrelhod. Hg. 7, diovre 
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Fig. 7. T he localizaSion of initial ncr tc ur for th 



the result of the determination of the iniha] snake (black 
circle) superimposed on the onginal color retinal image of 
Fig. 1. 

Because of the snake locks onto homogeneous regions 
enclosed by strong gradient informahon the ROI pari: of the 
enhanced image shown in Fig 5 (b) is processed bj the 
"closing" moqihological operation (dilation then erosion) to 
remove the regions of blood vessels that are found on the 
optic disk to create fairiy constant region before applymg the 
snake method. Fig. 8, shows the result (N B The ROI pari 
is zoomed to show deariy the fairiy constant region after 
applying the morphological closing operation). 

Finally the active contour algorithm is applied for 
accurately locating the optic disk and the result (the optic disk 
contour) is shown in Fig. 9(b). Comparing this result with 
Fig. 9 (a) (which shows the "ROI" part of the original color 
retinal image shown in Fig. 1), it is noticeably that, a perfect 
contour is obtained that accurately well define the optic disk. 

IV. Conclusion 

TIk identification of the optic didi contour is an inpoitant 
task in many opJithalmic ^::(Jicalions. The color retinal image 
is fully automated till readiing the exact contour of the optic 
di^ The badiground is dininated, Xh&i the intenaty image 
is eihanced to increase its contrast. The location of the of)tic 
di^ oaiter is defined to locate a suitaHe region of intaest 




"ROI". The rest of the piocessing steps are apfJied only to 
this region. This reduces the fxocesang time, prevails from 
defining a false initial contour, and makes it suitaUe for the 
online ^plication^ e.g. retinal identification. In tiie "ROI" an 
initial dtcular aiake is automatically estimated to be used as 
an initial contour in the active contour rrelhod 
Morpdiological ope:ations are crucial b^ore ^::fJying the 
active contour method. Finally the exact contour of the optic 
didt is detamined 
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